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Pritchard’s  Approximation  in  Array  Modeling 


by 

C.  L.  Scandrett  (Department  of  Mathematics) 
S.  R.  Baker  (Department  of  Physics) 
Naval  Postgraduate  School 
Monterey,  CA 

Abstract 


An  investigation  into  the  applicability  and  accuracy  of  Pritchard’s  approximation  for 
closely  packed  transducer  arrays  is  undertaken.  A  new,  ’’modal”  Pritchard  approximation 
is  developed,  based  upon  normal  modes  of  the  acoustic  medium,  and  is  tested  for  arrays  of 
acoustically  hard  spheres  to  ascertain  its  accuracy  in  determining  the  mutual  acoustic  radi¬ 
ation  impedance  between  array  elements.  For  ka  «  1,  it  is  found  that  the  modal  Pritchard 
approximation  works  quite  well  in  approximating  the  mutual  radiation  impedance  of  a  two 
element  array,  even  for  relatively  close  spacing,  but  for  arrays  of  three  or  more  scatterers  in 
close  proximity,  the  approximation  may  have  relatively  large  errors.  The  effect  of  neglect¬ 
ing  inter-element  scattering  is  analyzed  for  the  monopole  to  monopole  scattering  of  various 
configurations  of  a  three  element  array  and  a  sixteen  element  double  fine  array. 
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1  Theoretical  Background 


The  so-called  Pritchard  [9]  and  “generalized  Pritchard”  methods  [2]  are  popular  tech¬ 
niques  used  to  approximate  the  acoustic  interactions  between  elements  in  an  underwater 
acoustic  array.  Origins  of  the  method  can  be  traced  back  to  an  asymptotic  analysis  for 
small  spherical  scatterers  in  an  acoustic  medium  performed  by  Karnovsky  [7],  who  in  turn 
references  (as  does  Pritchard  in  his  paper)  earlier  work  done  (including  baffled  pistons)  by 
Rayleigh  in  1903  [10]. 

In  Pritchard’s  paper,  the  approximation  is  based  upon  an  analytical  treatment  of  acoustic 
interactions  in  an  array  of  baffled  circular  pistons  loaded  by  a  semi-infinite  fluid  medium. 
Extensions  to  more  general  arrays  are  possible  however  if  the  following  three  criteria  are 
met: 

•  The  mutual  impedance  between  any  two  elements  of  the  array  can  be  considered  with¬ 
out  regard  to  other  array  elements  in  the  construction  of  a  global  impedance  matrix 
for  the  entire  array. 

•  The  dimensions  of  an  array  element  (in  our  case,  the  sphere  radius  a)  is  small  relative 
to  the  acoustic  wavelength  (ka  <C  1,  where  k  =  ui/c  is  the  wavenumber  of  the  time 
harmonic  problem  and  c  is  the  acoustic  wave  speed). 


The  ratio  of  array  element  dimension  to  inter-element  distances  (d)  is  small  (—  <§;  1). 

(JL 


While  it  is  true  that  if  the  second  and  third  conditions  are  met,  so  must  the  first,  the 
method  is  often  applied  when  the  latter  two  conditions  are  not  strictly  satisfied,  with  very 
good  results.  This  is  most  likely  due  to  the  first  condition  being  met.  In  Pritchard’s  paper  the 
first  of  these  approximations  is  exactly  satisfied,  since  his  array  elements  are  coplanar  with 
the  baffle.  When  the  assigned  velocity  of  all  but  one  of  the  pistons  vanishes,  the  remainder 
of  the  array  becomes  essentially  part  of  the  baffle,  and  the  exact  mutual  radiation  impedance 
between  it  and  any  of  the  remaining  array  elements  can  be  found  analytically  without  regard 
to  the  other  array  elements,  because  there  is  no  scattering. 

The  Pritchard  approximation  for  the  mutual  radiation  impedance  between  two  array 
elements  (for  example  pistons  1  and  2)  can  be  formulated  mathematically  as  follows.  A 
series  representation  of  the  mutual  acoustic  radiation  impedance  between  the  two  pistons  to 
leading  order  for  small  ka  and  small  a/d  is: 


_  _  rsinfcd  .cosfcdl  _ 

Z\ 2  ~  R22  b  *  —  =  -^22^0  (kd) 

in  which  Z12  is  the  mutual  radiation  impedance  on  piston  1  due  to  a  radiated  pressure 
from  piston  2,  is  the  zeroth  order  spherical  Hankel  function  of  the  second  kind  (an  e**1* 
time  dependence  is  used  throughout),  and  R22  is  the  real  (resistive)  part  of  the  self 
radiation  impedance  of  piston  2.  Pritchard  demonstrated  that  the  above  approximation 
is  quite  accurate  in  comparisons  with  exact  calculations  over  a  range  of  center  to  center 
distances  between  pistons  in  the  baffled  array. 
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Benthien  [2]  reported  on  extensions  of  Pritchard’s  “method”  to  other  types  of  arrays.  In 
particular,  Benthien  found  excellent  agreement  in  comparing  a  Pritchard  approximation  to 
experimental  as  well  as  numerically  determined  results  based  upon  the  boundary  element 
technique  CHIEF[13],  for  a  three  element  array  of  flextensional  transducers.  In  Benthien’s 
methodology,  the  mutual  radiation  impedances  between  array  elements  are  found  by  numer¬ 
ically  determining  the  free  field  self  radiation  impedance  of  a  single  array  element  ( Zseif ). 
This  uses  the  code  CHIEF.  The  resulting  matrix  is  then  used  to  obtain  the  mutual  radiation 
impedance  of  that  array  element  on  a  second  (identical)  array  element  by  multiplying  the 
self  radiation  resistance  matrix  ( Rseif )  by  h^\kd)  where  d  is  the  distance  between  the  two 
array  elements. 

A  mere  complicated  and  perhaps  more  accurate  calculation  would  apply  a  Pritchard  type 
approximation  on  a  finite  element  by  finite  element  basis  rather  than  from  array  element  to 
array  element ,  but  is  not  done  for  two  reasons.  One  reason  is  that  the  calculations  involved 
in  finding  the  mutual  radiation  impedance  between  array  elements  for  a  large  array  would 
be  greatly  increased.  A  second  reason  is  based  upon  the  underlying  assumption  of  the 
Pritchard  approximation  itself,  namely  that  a<d  where  a  is  a  characteristic  length  of  the 
array  element.  This  implies  that  finite  element  to  finite  element  spacing  between  two  array 
elements  whose  center  to  center  spacing  is  d  should  be  very  nearly  d  as  d  ±  o  «  d. 

In  the  Benthien  formulation,  the  mutual  radiation  impedance  matrix  so  constructed  is 
applied  to  the  velocity  degrees  of  freedom  (one  for  each  surface  finite  element  facet)  for  each 
element  of  the  array  to  produce  a  pressure  forcing  function  on  each  of  the  array  element 
fluid  loaded  surfaces.  One  can  then  solve  for  the  unknown  values  of  the  velocities  by  matrix 
inversion  techniques.  Following  Benthien,  one  could  also  normalize  the  matrix  equations 
utilizing  the  in  vacuo  eigenvectors  of  the  mass/stiffness  finite  element  matrices  of  the  array 
element  prior  to  solving  the  system. 

The  work  of  Blottman  et  al  [3]  should  be  mentioned.  These  authors  determined  the  mu¬ 
tual  radiation  impedance  between  a  two  element  array  of  transducers  numerically  wherein 
the  in  vacuo  normal  modes  of  the  transducers  were  used  as  degrees  of  freedom  of  the  sys¬ 
tem.  In  this  regard,  Blottman’s  work  resembles  that  of  Benthien’s  “normalization” .  Unlike 
the  work  of  Benthien,  it  is  assumed  that  inter-mode  (cross-)  coupling  of  the  self  radiation 
impedance  for  a  single  transducer  is  negligible  (i.e.  in  vacuo  normal  modes  don’t  couple  when 
heavy  fluid  loading  is  applied).  Modal  mutual  impedances  between  array  elements  are  found 
by  employing  a  boundary  element  technique  (EQI/ATILA)[4j  which  (unlike  a  Pritchard  ap¬ 
proximation)  numerically  accounts  for  scattering  effects  to  all  orders.  While  the  work  by 
Blottman  et  al  would  appear  to  have  little  to  do  with  the  Pritchard  approximation,  the  idea 
of  applying  a  “modal”  mutual  radiation  impedance  is  of  interest.  In  particular,  this  report 
considers  the  possibility  of  a  modal  Pritchard  approximation  which  might  improve  upon  the 
Pritchard  approximations  currently  in  use. 
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2  The  mutual  radiation  impedance  matrix  and  the 

Pritchard  approximation 

Given  a  fluid  loaded  array  of  M  elements/transducers,  and  that  the  total  pressure  on 
array  element  n  is  represented  by  Pn,  its  value  can  be  formally  written  in  terms  of  the  mutual 
and  self  radiation  impedances,  and  the  prescribed  normal  surface  velocities  on  each  of  the 
array  elements  (Vm  for  array  element  m),  by  the  sum 

M 

Pn  ~  ^  ^  Znm  Vm 

m=  1 

These  pressures  can  in  turn  be  used  in  conjunction  with  surface  velocities  to  calculate  pres¬ 
sures  at  arbitrary  field  points  in  the  acoustic  medium. 

Crucial  to  Pritchard’s  method  is  an  accurate  representation  of  the  self  radiation  impedance 
of  an  array  element  since  it  is  the  basis  for  determination  of  all  mutual  radiation  impedances 
imposed  by  that  element  on  other  elements  in  the  array  (it  is  used  to  find  the  Znm  along 
that  row  of  the  global  impedance  matrix).  If  the  self  radiation  impedance  is  approximated 
by  the  free  field  radiation  impedance  for  the  array  element,  the  effect  of  inter-element  array 
backscattering  is  necessarily  neglected.  Backscatter  might  however,  be  taken  into  account 
at  some  level.  At  one  extreme,  one  could  numerically  find  the  self  radiation  impedance  of  a 
single  element  which  accounts  for  all  scattering  between  array  elements,  but  that  would  be 
tantamount  to  solving  the  full  array  problem  numerically,  precluding  the  need  for  a  Pritchard 
approximation  altogether.  (Incidentally,  this  is  exactly  what  Pritchard  was  able  to  do  for  his 
baffled  array.)  A  second  alternative,  neglects  all  backscattering  except  what  is  between  those 
two  array  elements  for  which  the  mutual  radiation  impedance  is  being  found.  For  a  two  body 
array,  this  and  the  former  methods  of  finding  the  self  radiation  impedance  are  equivalent.  A 
difficulty  with  this  latter  technique  is  that  the  self  radiation  impedance  expression  used  in  a 
Pritchard  approximation  for  each  array  element  would  depend  upon  the  geometry  (distance 
and  orientation)  between  it  and  each  member  of  the  array.  Calculations  would  then  have  to 
be  done  on  every  combination  of  two  element  arrays  possible  in  the  full  array,  leading  to  a 
considerable  amount  of  work  for  “large”  arrays. 

If  the  free  field  self  radiation  impedance  is  used  in  the  Pritchard  approximation,  and 
the  velocity  on  the  surface  of  each  array  element  is  known,  the  field  pressure  at  any  point 
in  the  acoustic  medium  can  be  found  by  a  simple  superposition  of  the  radiated  pressures 
resulting  from  each  of  the  array  elements.  The  surface  pressure  on  the  nth  array  element 
can  be  expressed  as  the  sum  of  its  self  radiated  pressure  (ZnnVn)  plus  contributions  of 
radiated  pressures  from  all  other  array  elements,  through  the  mutual  radiation  impedance 
matrix.  Pritchard’s  approximation  essentially  substitutes  the  contribution  to  the  radiation 
impedance  from  the  mth  array  element  on  the  nth  element  by  a  point  source  with  amplitude 
equal  to  the  resistive  part  of  the  mth  array  element’s  self-radiation  impedance.  This  can  be 
written  as 
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M  M 

Pn  ~  'y  ]  Znm  Fm  ~  ^rmfree^n  "b  )  ]  Prrvmh 0  (kdnm)  Vm 

m=  1  m=l,m^n 

where  M  is  the  number  of  elements  in  the  array,  and  dnm  is  the  distance  between  array 
elements  n  and  m. 

To  gauge  the  accuracy  of  the  Pritchard  approximation,  an  exact  lorm  of  the  mutual 
radiation  impedance  between  elements  of  an  array  is  needed  which  is  more  general  than  the 
“flat”  array  considered  by  Pritchard  in  his  original  work.  An  analytical  form  for  the  mutual 
radiation  impedance  between  two  acoustically  hard  spheres  is  given  by  New  and  Eisler  [5] 
which  employs  a  Green’s  function  analysis  coupled  with  spherical  addition  formulas  [8].  Self 
and  mutual  acoustic  radiation  impedances  between  the  spheres  are  found  and  represented 
by  infinite  series  involving  spherical  Bessel  functions  and  spherical  harmonics.  Results  are 
then  given  over  a  range  of  center  to  center  distances  for  cases  in  which  one  or  both  of  the 
spheres  are  oscillating  in  a  monopole  fashion. 

The  surface  radiation  impedance  used  by  New  and  Eisler  is  based  upon  the  definition 
[14] 

Zrj  =  / f  P(r7>*(rj)dSj 

where  p(r )  and  v(r)  are,  respectively,  the  pressure  and  normal  velocity  on  the  surface  (sj), 
and  where  v(r)  is  further  defined  as  v(r)  =  V /3(f)  in  which  V  is  an  amplitude  applied  to  the 
spatially  varying  function  /3(f) .  (Asterisks  refer  to  complex  conjugates.) 

It  is  a  simple  matter  to  find  values  for  the  self  and  mutual  radiation  impedance  for  the 
two  element  array  analyzed  in  New  and  Eisler’s  work  since  the  two  spheres  are  identical.  For 
this  array,  the  total  radiation  impedance  seen  by  sphere  1  is: 


The  scaled  radiation  impedance  on  sphere  one  when  both  spheres  have  identical  velocities,  is 
given  by  the  sum  Zn+Z^,  while  Zn  is  found  by  calculating  the  surface  radiation  impedance 
of  sphere  one  while  assigning  a  zero  velocity  on  sphere  two.  To  find  the  mutual  radiation 
impedance  of  sphere  two  on  sphere  one  (Z12),  the  difference  between  these  two  radiation 
impedances  is  taken  (ZTl  \Vl=V2=1  —  Zri  |Kl_i  V2=0).  (Equivalently,  Z12  could  be  found  by 
determining  the  pressure  on  sphere  1  when  V\  =  0  and  V2  =  1.) 

In  reporting  their  results,  New  and  Eisler  make  the  following  observation  regarding  the 
two  sphere  array  and  the  radiation  impedance  on  sphere  1  ( Zri ): 

. . .  the  contribution  of  the  mutual  term  Z12  is  certainly  significant.  However, 
the  primary  contribution  to  the  deviation  of  Zri  from  the  free  field  value  is  not 
scattering,  but  merely  a  consequence  of  the  added  pressure  field  from  sphere 
number  2.  One  can  conclude,  therefore,  that,  in  the  transition  region  (ka\  ~ 
ka,2  ~  1),  the  effects  of  scattering  on  the  radiation  impedance  of  one  spherical 
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transducer  in  the  presence  of  a  second  spherical  transducer  are  barely  significant, 
even  for  close  spacings. 

This  would  appear  to  validate  a  “Pritchard”  type  approximation  for  an  array  of  spheres,  up 
to  and  including  fe«l,  and  small  inter-array  element  spacings. 

In  our  analysis,  spherical  harmonic  functions  and  the  spherical  addition  formula  are  used 
to  find  an  appropriate  “modal”  Pritchard  approximation  based  upon  normal  modes  of  the 
acoustic  medium.  Our  work  explores  the  nature  of  this  approximation  with  comparisons 
to  exact  solutions  found  by  a  T-matrix  formalism  [11],  in  which  scattering  effects  between 
array  elements  are  fully  taken  into  account. 

3  Derivation  of  the  Modal  Pritchard  Approximation 

The  time  harmonic  radiated  or  scattered  pressure  field  is  represented  in  terms  of  out¬ 
going  spherical  Hankel  functions  applied  to  the  set  of  spherical  harmonics.  For  a  single 
radiator/scatterer  with  a  local  coordinate  system  written  with  the  index  “j”,  the  functional 
form  of  the  pressure  field  (jpj)  is  written 

N  n 

Pj{rji  fij)  =  X^  X/  ^rnn^ri(krj)Qn  ( 9j ,  <f>j) 

n=0  m—~n 

where 


-  K(<xx<>Ym* 

are  spherical  harmonic  functions. 

Representation  of  a  single  outgoing  spherical  wave  using  the  spherical  addition  formula, 
is  given  by  [8] 


OO  V  n+l/  1 

Mfe-»)STO>.&)  =  £  £  £ 

p—  |n— 1/| 
p>\m-ix\ 


a(u,p,  n,  p,  m)ju(kr1<)hp(kri>)n£+m(612 ,  0i2)fWi>  <j> i) 


(the  prime  on  the  summation  over  p  indicates  jumps  of  2  in  the  sum) 
(rj,  0i,  <f>i )  =  spherical  coordinates  relative  to  system  1 
(t*2,  #2,  <f> 2)  =  spherical  coordinates  relative  to  system  2 
(ri2>  #12, 0i2 )  =  origin  of  system  2  relative  to  1 
n>  =  max{r2,r12}  ru  =  min{r2,ri2} 
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and  where  the  coefficients  a(. . .)  are  related  to  the  Wigner  3-j  symbols  used  in  quantum 
mechanics. 

Derivation  of  the  modal  Pritchard  approximation  begins  with  a  description  of  a  two 
element  array.  Spheres  have  common  radius  a,  and  velocities  specified  on  spheres  1  and  2 
are  given  by  individual  normal  modes: 

Vt  =  —  ,*) 

up  dr  ri=a 

v*  =  ~z%  = 

r 

Letting 


and  noting  that 

SCOT  -  feTT  -  «  =  (-i)”+mn” 

the  pressure  from  sphere  2  translated  to  sphere  1  coordinates  (with  r\  <  r  12)  is  given  by: 


CO  *2 


£2=0  S2=— *2 


P2  =  E  E  /*s2t2  {  E  E  3vi{kri)  E '  Pi  1  S2)h<7,(kd)Q^l(e1, <j>i)Q,s0 


1/1=0  M=-Vi 


<Tl=\t2-V\\ 

CT1>|S2-Mll 


Similarly,  the  pressure  from  sphere  1  translated  to  sphere  2  coordinates  is 


OO  t\  f  OO  1/2 

pi  =  E  E  i  E  E  >2(^2) 

*1=0  51— -*1  (l/2=0  Ai2  =  -^2 


tl+1^2 

<72=|tl-l^| 

<^2>kl-/i2| 

Upon  application  of  the  boundary  conditions: 


KJ’fCff,,*) 

K£!tC(02,<« 


and  formally  applying  normal  modes  with  division  by  derivatives  of  Hankel  functions  applied 
at  ka  leads  to  the  system  of  equations 
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-ipcV™1 

h'ni(ka) 


$mmi  & nn  i 


&nn-L 


+ 


ini  (^Q)  V'  V'  ^2 

K  Wii(2  S2(2 


^2+^1 

E  '  a(ni,  <ri,  t2,  mh  s2)hffl  {kd)^~mi 

<Ti=\t2-ni\ 

-<’,1>|S2— nii| 


~jpcV£2 

h'n2{ka) 


^mm2^nn2 


—  ■^■mn^mm2^nn2 


+ 


in2(fc°) 

hn2(ka) 


OO  ti 

E 


tl=0 


E  '  a(n2,  (t2,  tu  rn2,s1)ha2(kd)QsJ2-m2(-l)a2+Sl-m2 

^2=\ti-n2\ 

,<T2>|si— Tn2| 


(The  terms  above  are  Kronecker  delta  functions  which  are  one  when  i  =  j  and  zero 
otherwise.) 

Truncating  the  series,  letting 


Bln 

77171 


-ipcV^1  c 

h>ni(ka)  *""* 


Jnn\ 


R2 

°  mn 


~ipcV£2  c  , 

h;2(ifca)  mm2<5nnz 


and  rewriting  these  equations  in  matrix  form  (e.g.  A1  =  [Aq0,  Ain,  A[0,  A\x,  A]_22,  •  •  -]r) 


I 

k12  ' 

' A 1  ' 

B1 

,k21 

I 

A2  _ 

B2 

where  the  [(si,  ti),  {s2,t2)}  entry  of  the  K\2  matrix  is 


■Kl2((«l,*l),(S2,*2)) 


j>  (ka)  t2+tl 

T777-T  IT'  o(<i)o-i5i2,s1,s2)/i<Tl(A:d)fi*2-Sl 

<Tl>|52—Sl| 


and  the  [(s2,t2),  (si,ti)]  entry  of  the  K2i  is 


K21  {(s2,t2),  (si,ti)) 


fs2(ka) 

h'S2  (ka) 


ti+t2 

E  '  ®(^2>  <72,  <1,  52,  Sl)^2(fcd)fi^-S2(-l)^+^-S2 


(72=|tl-t2| 

<^2>|si-S2| 


Solution  of  the  above  matrix  equation  for  the  unknown  amplitudes  provides  an  “exact” 
answer  to  the  array  problem,  and  is  equivalent  to  the  solution  found  by  employing  the 
T-matrix  formalism  of  Scandrett  and  Baker  [12]. 

For  fca<l  asymptotic  forms  for  the  ratio  of  the  spherical  Bessel  functions  are  [1] 
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j^jko)  (  0[(fca)3]  for  n  =  0  or  1 

h'n(ka)  ~  \  0[(ka)2n+1]  for  n  >  2 


as  ka  — >  0 


Therefore  the  magnitude  of  the  matrix  entries  Ai2  and  A2i  are  of  order  (ka)3,  and  an 
approximate  inverse  to  the  matrix  can  be  found  which  is  accurate  to  order  (ka)6 


+  0[(ka)6] 

and  an  explicit  calculation  for  the  unknown  scattering  amplitudes  is  possible: 

A1  «  B1  —  Kn B2  and  A2  «  B2  -  A21B1 

To  obtain  the  radiation  impedance  component  on  the  (ni,mi)  mode  of  sphere  1  given 
that  sphere  2  is  radiating  in  the  (n2,  ra2)  mode,  the  pressure  amplitudes  are  substituted 
into  the  radiation  impedance  integral  formula,  and  all  terms  of  order  (ka)6  or  higher  are 
neglected.  One  is  left  with  the  expression: 


I 

K12' 

p  O' 

+ 

0 

K12' 

r1- 

1 

-Kn' 

_K2i 

I 

"1 

°  I 

_k21 

0 

/  - 

.  -Kn 

I 

IL 


7  _  ~1pc 

“  l/m, 

.  ni 

_ c  3n  1  (ka)  -wrrni 

ni°h'ni(ka) 

a>|m2— mi| 


|Mlft0(^]2  +  3cos(2^)|M[ftl(M)]2 


+  0[(*g)6] 


where  H (4  —  —  n2)  is  a  Heaviside  function  which  is  zero  for  n\  +  n2  >  4  and  equals  one 

otherwise. 

The  middle  term  in  the  above  expression  is  the  sole  contribution  to  scattering  from  one 
sphere  to  the  other,  and  is  present  only  when  one  of  the  spheres  is  radiating  with  a  nonzero 
amplitude  in  the  monopole  mode  (at  least  to  0[(/ca)6]).  The  two  parts  of  this  term  result 
from  monopole  and  dipole  backscatter  from  sphere  two  which  affect  the  amplitude  of  the 
scattered  monpole  pressure  amplitude  on  sphere  one.  The  backscatter  terms  are  to  leading 
order  (ka)5,  and  if  neglected,  one  obtains  an  approximation  for  the  radiation  impedance  to 
0[(&a)5]  which  completely  neglects  all  scattering  effects,  and  is  consistent  with  using  the 
free  field  radiation  impedance  for  the  self  radiation  impedance,  as  outlined  in  the  previous 
section. 
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JT  \ 


— %pc 
Ymi 

V7li 


-iff  (3  -n!-n2) 
(ka)2h'ni{ka)h'n2{ka)  ”2 


n2+ni 

'  a(n1,a,n2,m1,m2)^(W)fl”rrai  > 

<7=1712-0!  | 

<7>|77l2— 77ll|  , 


The  first  term  which  is  multiplied  by  V™1  is  the  self  radiation  impedance 


and  equals  the  free  field  radiation  impedance  of  sphere  one.  The  summation  term  multiplying 
V™2  is  the  mutual  radiation  impedance  of  sphere  2  on  sphere  1,  and  can  be  thought  of  as 
an  approximation  to  the  translation  of  the  radiated  pressure  from  sphere  2  onto  sphere  1: 


^12(n2,m2)— (ni,mi)  &  JJ$  H^i1  (^1>  ^l)l|2  ^S1 

—iff  (3  —  ni  —  712)  _ 

(ka^h^ika^ika)  mi 

0*>|rrt2 — T7li  | 

The  matrix  form  for  the  non-neglected  terms  in  the  modal  Pritchard  approximation  to 
the  mutual  radiation  impedance  is: 


Z(00)  (00) 

Z(00)(l-1) 

%>0)(10) 

^(00)(11) 

^(00)(2-2) 

-Z(00)(2-l) 

Z(00)  (20) 

^(00)(21) 

^(00)(22) 

0  ...  “ 

^(1-1)(00) 

z(i-i)(n) 

0 

0 

0 

0 

0 

0  ... 

Z(10)(00) 

Z(10)(l-1) 

Z(10)(10) 

Z(10)(ll) 

0 

0 

0 

0 

0 

0  ... 

Z(H)(00) 

z(ll)(l-l) 

Z(U)(10) 

z(ii)(ii) 

0 

0 

0 

0 

0 

0  ... 

<Z(2-2)(00) 

0 

0 

0 

0 

0 

0 

0 

0 

0  ... 

^(2-l)(00) 

0 

0 

0 

0 

0 

0 

0 

0 

0  ... 

Z(20)(00) 

0 

0 

0 

0 

0 

0 

0 

0 

0  ... 

2’(21)(00) 

0 

0 

0 

0 

0 

0 

0 

0 

0  ... 

^(22)(00) 

0 

0 

0 

0 

0 

0 

0 

0 

0  ... 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0  ... 

It  is  interesting  to  compare  the  real  part  of  the  self  impedance  to  the  amplitude  of  the 
mutual  impedance  in  the  special  case  that  both  spheres  are  radiating  identically  (B1  =  B2 
or  ni  =  n2  =  n  and  =  m2  =  m).  In  this  case,  the  real  part  of  Zn  =  Rn  becomes 


^  -  peg  {£g} 

-  -peg  iig.^>1)i^^),||;;i(fa)||2 

f ka)2n+ 2  rr 

[(n  +  1)(1)(3)  •  •  •  (2n  -  1  )^PC I JSl 
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while  Z\2  has  the  asymptotic  form 


[(n  + 


DdiTo-i WpcIL 


x  H(3  —  2n)^2'  a(n,a,n,m,m)ho{kd)£ft0 

<7=0 


When  only  the  first  term  in  the  series  for  Z12  (the  a  —  0  case  which  results  in  an  ho(kd) 
term)  is  kept,  the  “standard”  Pritchard  approximation  results.  In  this  instance  note  that 
a(n,  0,  n,  m,  m)  =  1  for  all  n  >  0  and  —n<m<n.  The  remaining  terms  in  the  series  have 
an  angular  dependence  which  affects  the  amplitude  of  Z\i  as  the  relative  angles  between 
array  elements  vary.  Such  an  angular  dependence  is  not  built  into  the  standard  Pritchard 
approximation  -  only  the  relative  distance  between  elements  as  expressed  by  the  factor  ho(kd) 
is  present. 

4  Results 


4.1  Two  body  interactions 


In  this  section  several  graphs  will  be  displayed  in  an  effort  to  demonstrate  the  modal 
Pritchard  approximation  given  above.  The  first  set  of  numerical  experiments  involve  the 
New  and  Eisler  problem  of  an  array  of  two,  identical,  acoustically  hard  spheres  separated  by 
a  distance  d.  Each  sphere  has  radius  a  =  1/2,  and  is  assigned  a  modal  velocity  in  such  a 
way  that  the  mutual  radiation  impedance  can  be  found. 

Besides  the  modal  Pritchard  method,  there  are  two  other  approximations  displayed.  The 
results  labelled  “Pritchard’s  method”  and  “Simplified  Pritchard”  are  found  as  follows.  For 
the  “Pritchard  method”,  an  “exact”  value  for  R\\  which  includes  all  backscattering  is  used 
in  the  approximation,  while  the  “simplified  Pritchard”  employs  only  the  free-field  values  of 
the  self  radiation  resistance  term  Ru,  and  therefore  neglects  backscatter. 

The  first  series  of  graphs  display  the  monopole  to  monopole,  axially  symmetric  dipole 
to  dipole,  and  axially  symmetric  quadrupole  to  quadrupole  mutual  radiation  impedances 
for  the  two  spheres  which  are  aligned  with  the  z  axis.  The  final  graph  of  quadrupole  to 
quadrupole  is  included  only  to  show  the  accuracy  of  the  modal  Pritchard  approximation 
even  through  the  quadrupole  term,  in  spite  of  the  fact  that  this  term  is  of  order  ( ka )6  in  the 
expansion  for  small  ka  of  the  mutual  radiation  impedance  and  therefore  would  be  neglected 
in  an  application  of  the  derived  modal  Pritchard  approximation. 
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Z12  for  hard  spheres,  ka=1,  monopoles,  scaled  by  4rz  a2pc 


Z12  for  hard  spheres,  ka=1 ,  dipoles,  scaled  by  4 n  a2pc 


I  i  i  i  i  i  i - 1 - 1 - 1 - r 
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It  would  appear  from  the  above  graph  that  backscattering  effects  in  the  i?n  determi¬ 
nation  are  negligible  by  comparing  results  of  the  “simplified  Pritchard”  and  “Pritchard” 
methods.  Additionally,  the  “modal  Pritchard”  seems  to  very  closely  match  the  “exact” 
spherical  addition  formula  results  which  includes  all  scattering  effects.  The  results  of  the 
monopole  to  monopole  mutual  radiation  impedance  highlight  what  appears  to  be  a  phase 
shift  between  the  “modal”  and  “simplified”  Pritchard  methods.  Because  of  their  analytical 
form,  a  comparison  between  Pritchard  and  the  exact  solution  can  be  analytically  performed 
by  comparing  the  simplified  Pritchard  with  the  modal  Pritchard.  The  formulas  for  the  modal 
and  simplified  Pritchards  (scaled  by  Aira?pc)  as  a  function  of  the  distance  between  spheres 
(i d )  is 


modal  Pritchard  monopole  to  monopole  = 
simplified  Pritchard  monopole  to  monopole  = 


ho(kd ) 


-1 


hQ{kd)lZe 


[kah'0(ka)}2 
iho(ka) 


h'0(ka ) 


Taking  the  ratio  of  these  two  expressions,  one  has 


h0(kd ) 


-1 


kah'0(ka)Y 


ho(kd)lZe 


— iho(ka ) 


[  h'0{ka)  J 


fo(ka)  ±  Woi^  ; 

j'o(ka)  -  iy'o(ka) 


«  l  +  y(fca)3 


as  ka 


0 


One  can  also  see  from  the  above  expression,  that  the  amplitude  of  the  ratio  is  one  with 
a  constant  phase  shift.  When  ka  =  1,  this  phase  shift  is  about  24.6  degrees,  which  is  what 
is  seen  in  the  monopole  to  monopole  figure  above. 

For  the  dipole  to  dipole  case,  the  relative  position  of  the  spheres  introduces  an  angular 
dependence  in  the  mutual  radiation  impedance  predicted  by  the  modal  Pritchard  approxi¬ 
mation  which  is  not  present  in  the  simplified  Pritchard  approximation.  As  in  the  monopole 
to  monopole  case,  one  finds  that  the  simplified  Pritchard  approximation  is  essentially  equiv¬ 
alent  to  the  Pritchard  method  which  accounts  for  backscattering,  but  to  compare  with  the 
exact  and/or  modal  Pritchard  approximation,  the  results  must  be  magnified  by  a  factor  of 
three.  To  see  where  the  factor  of  three  comes  from,  one  can  again  compare  the  (scaled) 
modal  and  simplified  Pritchard  approximations 


modal  Pritchard  dipole  to  dipole 
simplified  Pritchard  dipole  to  dipole 


3  [kah[(ka)]‘ 


{h0(kd)-2h2(kd)Q°2(e12,<f>12)} 


ho(kd) 


He 


—ihi  ( ka ) 


L  hi(ka)  J 


If  one  assumes  kd  1,  the  approximation  h2(kd)  ~  —ho(kd)  can  be  made,  and  with  the 

fact  that  f2°(0i2)  =  +  3cos(2#i2)]/4,  the  ratio  of  the  two  expressions  can  be  written: 
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ho(kd)  — 


3[1  +  cos(2#i2)]/2 
[kah[(ka)]2 


ho(kd)Ke 


—ihi(ka) 
h[(ka ) 


3[1  +  cos(26>i2)]  j  j[{ka )  +  iy[(ka) ) 

2 


3[1  +  cos(2^i2)] 
2 


+  0[(fca)3]  as  ka 


0 


In  addition  to  the  phase  shift  embodied  by  the  ratio  of  the  spherical  Hankel  functions,  one 
sees  an  amplitude  variation  with  respect  to  the  angle  #i2  between  the  spheres  (in  this  instance, 
the  axis  of  the  dipole  aligns  with  the  z  axis  from  which  the  angle  812  is  measured).  For  the 
geometry  considered  in  the  graph,  the  angle  812  =  0  leading  to  a  factor  of  3  needed  to  bring 
the  Pritchard  approximation  in  line  with  the  analytic  solution.  The  amplitude  and  phase 
shift  are  unfortunately  not  the  only  differences  between  the  solutions.  For  small  distances 
d  the  approximation  used  for  h2(/cd)  is  no  longer  valid,  and  this  can  be  readily  seen  in  the 
graph  of  the  dipole  to  dipole  mutual  radiation  for  small  values  of  d/a. 

Below  is  the  axially  symmetric  quadrupole  to  quadrupole  mutual  radiation  impedances 
of  the  various  Pritchard  approximations  along  with  the  spherical  addition  result.  The  modal 
Pritchard  result  would  actually  be  zero  if  we  neglected  all  terms  of  order  ka5  since  quadrupole 
to  quadrupole  radiation  is  of  higher  order.  This  graph  is  only  included  to  show  the  continued 
success  of  the  approximation  embodied  by  the  formula  derived  in  the  previous  section  for  Z\2- 
Note  also,  that  the  magnitude  of  these  quadrupole  impedances  are  an  order  of  magnitude 
less  significant  than  the  dipole  terms. 
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4.2  Three  body  interactions 

One  might  hope  that  the  modal  Pritchard  method  works  well  even  for  situations  in  which 
there  are  more  than  two  elements  in  a  given  array.  This  would  indeed  be  true  if  ”  third  party” 
scattering  were  negligible.  Third  party  is  meant  to  imply  single  or  multiple  scattering  of  an 
incident  pressure  from  one  array  element  to  another  (thereby  affecting  the  mutual  radiation 
impedance  between  the  two),  by  way  of  scattering  from  a  third  obstacle  or  array  element. 
This  is  an  important  effect  when  considering  closely  packed  arrays  of  transducers.  If  one 
considers  only  two  element  arrays  in  the  production  of  mutual  radiation  impedances  between 
array  elements,  such  scattering  effects  are  necessarily  neglected. 

The  most  elementary,  and  perhaps  most  enlightening  case  to  consider,  are  three  element 
arrays,  in  which  one  of  the  elements  acts  solely  as  an  acoustically  hard  obstacle,  and  is  moved 
around  relative  to  the  remaining  two,  which  are  fixed.  For  such  an  array,  one  can  anticipate 
the  negligible  effects  of  the  third  array  element  when  it  is  far  from  the  two  active  elements, 
and  one  can  quantify  the  effect  it  has  on  the  mutual  radiation  impedance  when  it  is  moved  to 
the  nearfield.  Because  the  two  active  array  elements  are  fixed  in  space,  the  modal  Pritchard 
approximation,  (and  for  that  matter  all  of  the  Pritchard  approximations),  produce  constant 
mutual  radiation  impedances  between  the  two  elements,  regardless  of  the  positioning  of  the 
third.  Because  of  its  accuracy  in  the  two  element  array  problem,  only  the  mod5*1  Pritchard 
approximation  will  be  compared  with  the  spherical  addition  results,  to  see  h<  mportant 
third  party  scattering  is  in  the  determination  of  mutual  radiation  impedances. 

Z13  for  hard  spheres,  ka=1,  monopoles,  scaled  by  4na2pc 


scaled  distance  of  second  sphere  from  array  axis  D/a 
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In  the  first  test  case,  the  two  active  array  elements  (numbered  elements  1  and  3)  are 
aligned  to  the  z  axis,  and  are  at  a  constant  center  to  center  distance  of  A/2,  which  for  ka  =  1 
implies  a  separation  of  7r/2.  The  third  element  (element  number  2)  of  the  array  starts  at 
a  distance  of  A/4  away  from  the  axis  of  the  array,  halfway  between  the  two  active  array 
elements.  It  then  moves  axially  away  from  the  two  active  array  elements.  Above,  is  the 
comparison  of  the  monopole  to  monopole  mutual  radiation  impedance  of  sphere  3  on  sphere 
1. 

Notice  that  there  is  a  sizeable  difference  in  the  resistance  and  reactance  terms  of  the 
mutual  radiation  impedance  when  sphere  number  2  is  in  close  proximity  to  the  radiating 
spheres.  In  fact,  in  the  third  graph  of  the  above  series,  the  relative  error  of  the  magnitude 
of  the  mutual  radiation  impedance 


Percent  error  equals  100  x  l|Zl3l|m°dal  Pritchard  ll^llspherical  addition 

I  ^13 1 1  spherical  addition 

shows  relative  errors  as  great  as  15%. 

The  deviations  of  the  mutual  radiation  impedance  amplitudes  can  be  explained  by  a 
relatively  simple  argument.  Depending  upon  the  total  distance  from  sphere  3  to  2  to  1,  the 
scattered  pressure  is  either  in  phase  (when  Z13  is  at  a  maximum),  or  out  of  phase,  with  the 
direct  pressure  from  sphere  3  to  sphere  1. 


|0o(ka){ho(kd13)-[ho(Kd12)]2io'(ka)/ho'(ka)}|| 
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This  occurs  when  the  total  scattered  distance  from  sphere  3  to  2  to  1  is  an  integral  multiple 
of  the  acoustic  wavelength.  This  may  seem  counter-intuitive  since  the  scattering  sphere  is 
acoustically  hard,  but  can  be  explained  by  a  simple  argument  involving  monopole  radiators 
and  the  spherical  addition  formula.  Considering  only  the  monopole  terms,  the  sum  of  the 
directly  incident  and  once  scattered  monopole  pressures  is  given  by  the  formula 

Aj0(ka)  |fto(M13)  -  (fco(Mi2)]2} 

where  d\z  is  the  distance  between  spheres  1  and  3,  c?i2  is  the  distance  between  spheres  1  and 
2  (and  also  the  distance  between  spheres  2  and  3),  and  A  is  the  amplitude  of  the  monopole 
pressure  radiated  from  sphere  3.  The  value  of  di2  can  be  represented  in  terms  of  sphere  2’s 
distance  from  the  z  axis  ( D )  by  dn  =  (A/4)2  +  D2.  In  the  figure  above,  the  amplitude  of 

the  above  factor  for  A  =  1  is  produced,  and  to  a  large  degree  explains  the  variation  seen  in 
the  amplitude  of  the  mutual  radiation  impedance  term  (Zn)  found  in  the  previous  graph. 

A  second  three  element  array  is  also  considered,  in  which  the  active  elements  are  the 
same  as  before,  but  in  which  the  scattering  element  moves  in  an  elliptical  path  around  the 
two  active  elements  in  such  a  way  that  the  two  active  elements  are  foci  of  the  ellipse.  In  this 
case,  the  scattering  distance  from  sphere  3  to  2  to  1  has  the  constant  value  3 A/2.  In  terms  of 
the  preceeding  three  body  graph,  the  point  which  would  intersect  the  elliptical  path  occurs 
when  D/a  «  4.44,  or  when  the  relative  error  in  Ziz  is  at  its  maximum  value  (corresponding 
to  the  Pritchard  approximation  being  greater  than  the  spherical  addition  result). 


Z1 3  for  hard  spheres,  ka=1 ,  monopoles,  scaled  by  4rca2pc 
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Graphs  of  the  resistamce  and  reactance  as  well  as  the  relative  error  in  the  amplitude  of 
are  given  for  angles  0  to  180  degrees.  At  0  degrees,  the  center  values  for  spheres  1,  2,  and  3 
are  respectively  0,  A,  and  A/2,  while  at  180  degrees,  the  z  values  of  these  spheres  is  0,  —A/2, 
and  A/2. 

In  this  instance  the  variation  due  to  scattering  is  nearly  uniform,  with  maximum  devia¬ 
tions  occuring  at  the  poles  and  perpendicular  to  the  “active”  array  axis.  Unlike  the  previous 
three  body  case,  monopole  scattering  alone  does  not  explain  the  variation  of  the  amplitude  of 
Ziz-  The  graph  below  shows  the  result  of  limiting  the  number  of  terms  retained  in  the  spher¬ 
ical  addition  formulation  to  monopole  only,  monopole  and  dipole  only,  and  monopole,  dipole 
and  quadrupole  terms.  The  graph  illustrates  that  the  monopole  to  dipole  back  to  monopole 
scattering  is  the  primary  reason  for  the  deviations  seen  in  the  amplitude  of  the  Z\z  term 
(the  monopole-monopole-monopole  and  monopole-quadrupole-monopole  terms  are  nearly 
constant  throughout  the  range  of  angles).  Note  also  that  the  result  of  keeping  only  through 
the  quadrupole  terms  very  nearly  matches  the  result  of  keeping  harmonics  up  through  order 
6. 


Z13  for  hard  spheres,  ka=1 ,  monopoles,  scaled  by  4rca2pc 


4.3  Double  line  array  of  hard  spheres 

As  a  final  example,  consider  a  double  line  array  of  hard  spheres  which  are  spaced  A/2 
apart  along  the  lines,  and  which  has  the  spacing  A/4  between  the  two  lines.  The  double 
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line  of  spheres  is  parallel  with  the  z  axis.  A  graphic  of  the  array  is  shown  below,  with  the 
numbering  of  the  spheres  starting  at  the  bottom  left,  and  ending  with  the  16th  element  at 
the  top  right. 


A\  A., Z 

Til 2  m  I  f' 

yrt2"- 


wavelength 
X=n  m 


9 's  ■ 


Y 

— - > 

Sphere  radius 
0.5  meters 


Til 4  m 

Orientation  of  the  spheres 
in  the  16  element  array 


In  the  first  graph  below,  the  source  level  (in  dB)  of  the  array  given  that  each  sphere  is  oscil¬ 
lating  in  a  breathing  mode  with  surface  velocity  one  is  shown.  The  spherical  addition  result 
is  compared  to  the  modal  Pritchard  calculation,  and  as  can  be  seen,  the  modal  Pritchai  i 
approximation  does  quite  well,  with  a  maximum  error  of  about  2  dB  at  broadside.  The  far 
field  results  are  scaled  by  the  far- field  amplitude  of  a  single  oscillating  sphere.  A  second 
graph  shows  the  magnitude  of  the  far  field  pressure,  scaled  by  a  single  radiating  sphere. 


Source  level  comparison  for  16  element  array,  hard  spheres,  ka=1 
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Source  level  comparison  for  16  element  array,  hard  spheres,  ka=1 


It  is  possible  for  this  array  to  compare  the  self  and  mutual  radiation  impedances  for 
the  monopole  terms  without  too  much  difficulty.  In  the  next  two  pictures,  the  self-radiation 
impedance  of  each  of  the  array  elements  is  given,  using  both  the  modal  Pritchard  approxima¬ 
tion  and  the  full  spherical  addition  results.  Recall  that  the  modal  Pritchard  result  neglects 
third  party  scattering  and  backscattering.  In  the  first  graph  of  the  self  radiation  impedance 
one  can  see  that  there  is  nearly  a  uniform  15%  relative  error  in  all  but  the  spheres  found 
at  the  extremities  of  the  array.  The  reason  for  the  variation  from  the  modal  Pritchard  is 
because  of  backscatter  not  being  taken  into  account,  and  that  at  array  elements  1  and  16, 
the  backscatter  is  less  than  for  any  of  the  remaning  array  elements.  The  second  graph  shows 
on  an  element  by  element  basis,  the  total  mutual  radiation  on  each  sphere  due  to  radiation  of 
pressure  from  all  other  array  elements.  Essentially,  it  is  the  sum  of  the  off  diagonal  elements 
in  the  monopole  impedance  matrix  of  the  array  in  any  given  row,  e.g.  for  array  element  i, 
this  value  is  essentially 

16 

‘ mutual  V'  r7 

Prod  =  2_,  Zi0 

j= 

since  all  array  elements  have  identical  (monopole  excited)  velocities  of  one. 
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5  Conclusions  and  recommendations 

The  work  reported  would  indicate  that  improvements  can  be  made  to  existing  methodolgy 
in  array  modeling.  In  particular,  for  low  frequency  active  array  calculations  with  close  array 
spacing,  it  would  appear  that  the  spherical  addition  formula  should  be  used  to  incorporate 
inter-element  scattering.  It  would  appear  however,  that  the  number  of  harmonics  needed  is 
only  through  the  quadrupole  term  for  accuracy  up  to  a  level  of  (ka)6,  and  so  if  a  T-matrix 
algorithm  is  adopted,  a  total  of  81  matrix  elements  should  properly  characterize  the  array 
problem. 

Great  accuracy  has  been  achieved  by  using  a  Pritchard  approximation  to  model  the 
mutual  radiation  impedance  between  array  elements  for  several  of  the  reasons  delineated  in 
the  above  report.  Improvement  would  seem  to  be  possible,  however,  by  adopting  a  modal 
type  Pritchard  approximation.  There  are  still  limitations  to  the  methodology,  however, 
which  should  be  addressed  in  an  actual  application,  such  as  the  need  to  model  third  party 
scattering.  The  manner  in  which  such  a  modal  approximation  might  be  incorporated  into  a 
finite  element  code  is  as  follows. 

The  finite  element  code  must  be  used  to  determine  the  self  radiation  impedance  of  a  given 
array  element.  Almost  any  boundary  element  code  would  do,  for  example  EQI  or  CHIEF, 
but  it  must  couple  to  the  finite  element  code  of  the  array  element.  Given  nodal  values  of 
the  normal  surface  velocity  (Vj)(or  iwx  the  normal  surface  displacement),  and  the  free  field 
self  impedance  radiation  matrix,  the  nodal  pressure  vector  (Pj)  on  the  jth  array  element  is 
given  by 


(Pj)  nodal  —  Zjj  Vj 

The  nodal  pressure  values  must  be  converted  to  modal  pressure  amplitudes  which  can  be  used 
by  the  spherical  addition  formula  for  translations  to  other  array  elements.  First  introduce 
an  intermediate  matrix  (A)  of  spherical  Hankel/spherical  harmonic  functions  evaluated  at 
the  surface  nodes  of  the  jth  array  element  (rs,  0S,  and  4>s) 


As,(n,m)  =  h{y(krs)Q™(6s,  <f>s) 

The  size  of  this  matrix  is  M  x  K  where  M  is  the  number  of  surface  nodes,  and  K  is  the 
number  of  harmonics  (9  if  use  only  harmonics  through  the  quadrupole  term).  The  equation 
to  be  solved  is 


(P)nodes  =  A(P) 


modes 


The  least  squares  solution  to  this  problem  is  found  by  solving  the  normal  equations,  resulting 
in 


(P)modes  =  (AHA)-1AH(P)nodes 
c 
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where  AH  is  the  Hermitian  conjugate  of  A,  and  C  is  the  sought  after  translation  matrix  from 
nodal  pressures  to  modal  pressures. 

Now  with  the  modal  radiated  pressures,  the  addition  formula  can  be  used  to  determine 
the  nodal  pressures  resulting  from  the  modal  radiated  pressures  through  an  application  of  a 
truncated  spherical  addition  formula.  The  translation  from  array  element  j  to  array  element 
i’s  nodal  coordinates  rt,  9t,  and  <f>t  is  given  by 

(sr“  = 

(h,mi  h,m2 

fa  +h  _ 

x  '  a(li,a,  l2,mi,m2)h<T(kd)£l™2~rni 

<r=|*2-*l| 

mi| 

where  d  is  the  distance  between  the  ith  and  jth  array  element,  and  H(  )  is  the  Heaviside 
function.  Introducing  the  matrix  G  which  is  M  x  K  (assuming  the  ith  array  element  also 
has  M  surface  nodes,  it  is  given  by 

K  h+h 

^  H(3  —  l\  —  h)jh (krtjQ™2 (9t, (f)t)  y  a(li,o’,l2,Tni,m2)hcr(kd)Qa2  1 

hi^2  a=\l2~h  | 

cr>\m2—mi\ 

Combining  the  three  matrices  into  a  modal  approximation  for  the  mutual  radiation 
impedance  leads  to 


Zij  ~  GnCtjZjj 

which  corresponds  to  the  modal  Pritchard  approximation  used  in  this  report,  and  which 
should  be  compared  to  the  standard  Pritchard  approximation  given  by 

Zij  ~  TZ’G  ( Zjj )  ho(kdij) 
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